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Diagrammatic techniques to compute perturbatively the spectral properties of Euclidean random 
matrices (ERM) in the high-density regime are introduced and discussed in detail. Such techniques 
are developed in two alternative and very different formulations of the mathematical problem and 
are shown to give identical results up to second order in the perturbative expansion. One method, 
based on writing the so-called resolvent function as a Taylor series, allows to group the diagrams 
in a small number of topological classes, providing a simple way to determine the infrared (small 
momenta) behavior of the theory up to third order, which is of interest for the comparison with 
experiments. The other method, which reformulates the problem as a field theory, can instead be 
used to study the infrared behaviour at any perturbative order. 

PACS numbers: 61.43.Fs, 62.10.+S, 

I. INTRODUCTION 

Random matrices [l| are N x N matrices whose entries are random numbers drawn from a certain probability 
distribution. Their statistical spectral properties in the large N limit describe a wide range of physical phenomena: 
nuclear spectra [H, quantum chaos [Jh, localization in electronic systems diffusion in random graphs @, liquid 
dynamics @ and the glass transition Q, complex networks Q, superstrings @. Random matrices may be grouped 
in a few universality classes according to their statistical properties For most of these classes, the density of 
eigenvalues follows Wigner's semicircle law. It has thus become of interest to identify ensembles where the semicircle 
law is modified in a non-trivial way. One such ensemble results when the corresponding physical problem has a 
conserved quantity (e.g. momentum in case of propagating excitations, or number density in diffusion problems). 
Under such circumstances, the random matrix that best describes the problem is typically a Laplacian matrix [5f, 
which has the property 



0. (1) 



This encodes the property that a vector whose components are identical is an eigenvector with eigenvalue zero. 

A kind of random matrices of particular relevance in the study of off-lattice systems are the so-called Euclidean 
random matrices (ERM) |2lj. Place N particles in positions Xi, i = 1,2, ... ,N, belonging to some region of 
D-dimensional Euclidean space, of volume V. The positions are drawn randomly from some probability distribution 
function P({xi}). The entries of an ERM are a deterministic function of these random positions, M^- = /(x, — Xj). 
If a conservation law is relevant for the problem at hand, we will rather have a Laplacian ERM: 

N 

My = kj /( x j - x fc) - f(*i - Xj)- (2) 
fe=l 

Note that we never find the same particle label twice in the argument of the function /(xj — Xj), since the term 
/(xj — Xj) cancels. In a diagonal term, 5jj/(xj — Xfe), the fcth particle shall be called a medium particle, while the ith 
particle will be the chain particle. 

The function / in Eq. @ is quite general: only rotational invariance and the existence of the Fourier transform 
f(p) are assumed (p = ^/p • p). Furthermore, even if in this work / will be a scalar function, for some applications it 
should rather be a matrix-valued function. It must be so, for instance, to account for the vector nature (longitudinal 
or transversal) of vibrational dynamics [l8j . Most of our results extend as such to this more general case. 
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ERMs describe topologically disordered systems, at variance with problems were the N positions {x^} are placed 
on a crystalline lattice |29{ . We will be considering a extreme case, in which the N positions are placed with uniform 
probability on the volume V. The particle-number density, p = N/V will be held fixed while we take the large N 
limit. Note that there are two sources of statistical correlation among the entries of matrix ©, even if the positions 
{xi} are totally uncorrelated. First, it is a Laplacian matrix, recall Eq. (JTJ). Second, due to the triangular inequality 
of Euclidean geometry, the distances from two neighbouring particles to a third one are necessarily similar. 

Specific applications of ERMs include disordered d-wave superconductors (Tlj . disordered ma gnet ic semiconduc- 
tors [12! (very similar to a spin-glass model [HI), Instantaneous Normal Modes in liquids |5, Tlij. vibrations in 
glasses M, the gelation transition in polymers and vibrations in DN A |20| . ERMs have been studied 

analytically and numerically both in the low particle- number density regime [3, [22|, l23j) and for high densities [HI EH - 

Mm- 

In order to compute the basic spectral properties of ERM it turns out to be convenient to introduce the resolvent, 
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N 

G(p, z) = lim — >^ 



e ip(xj-x ; 



M 



(3) 



where the complex number z = X + vq has a tiny imaginary part rj and the overline stands for an average over the 
{xj}. If the ERM describes physical excitations (phonons, electrons, etc.) in topologically disordered systems, the 
resolvent Q corresponds to the single- particle Green function, or propagator, for such excitations. If the system is 
isotropic, the resolvent depends only on p. 

The density of eigenvalues 5(A), or density of states (DOS), is given by 

g(X) =— lim Im[G(p,A + iO+)]. (4) 

7T p— >oo 

This limiting behaviour is characteristic of topologically disordered systems. It does not hold for lattice systems. We 
note as well that the constraint ([1]) implies that a plane wave e lp Xi is an eigenvector of the matrix @ if p = 0: 

G{0,z) = - (5) 

z 

As we shall discuss below, the resolvent takes a very simple form in the high-density limit (it is actually the bare 
propagator of the theory) : 

Go(p,z) = Vv <P)=p[f(°)-f(p)}- (6) 

z - e(p) 

The physical interpretation is quite appealing [loj . The system behaves as an elastic continuum medium. In the large 
p limit, the plane waves e lp Xi become exact eigenvectors of the matrix @, with eigenvalues given by the dispersion 
relation e(p) = p[f(0) — /(p)]. In particular, for small p, e(p) = c 2 p 2 + C(p 4 ), where c is the speed of sound. This 
neat physical picture motivates the introduction of a high-density expansion. 
At large, but finite p, the resolvent can be written 

G(p,z) = 1 -. (7) 

z - e(p) - E(p, z) 

The self-energy S(p, z), which is introduced to encode all the information about the interactions (a standard practice 
in the Green function formalism), vanishes when p tends to 00. In our case, the interaction involved is that between the 
propagating excitations and the topological disorder. An important theoretical challenge is to compute the self-energy 
at finite densities p. In fact, in this case an eigenvector can be thought as a packet of plane waves (see section [TTT) . 
The width of such packet is related to the imaginary part of S. 

Some of us have argued that in the limit of small p, z the leading term at 1/ p 2 order has the form fl5l — fl^zj ] 



D-2 



lmE{p,z + iO + ) =Az— p 2 + 0{z—p A ,z u p 2 ), (8) 

(D is the space dimension, while A is an amplitude) . This has been disputed recently by an independent computation 
to order l/p 2 , claiming that the actual analytic structure of the self-energy for small z and p is (30| 

lmS(p,z + i0 + ) cxzt p 2 + ... ; (9) 
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analogous to what one finds in the Rayleigh theory of scattering and in lattice models where disordered spring 
constants mimick the effect of topological disorder |29| ■ 

By reconsidering in detail the perturbative expansion, in this work we show that the prefactor A in Eq. is 
actually null, due algebraic cancellations, and that this cancellations arise at all orders in the perturbative expansion 
in l/p. This is not related to any known symmetry of the problem, but rather reflects the mathematical structure 
of the perturbative contributions. On the other hand, we will also show that the result in Ref. H3, recall Eq. ([9]), is 
incomplete, since the imaginary part admits a formal expansion for small z 



The constraint ([5]) implies that g n (0) = for all n, so that in general g n (p) = A n p 2 + 0(p A ). However, we find 
that, for all functions / and all p, Ao = 0, so that go(p) ~ p 4 while gi{p) ~ p 2 ■ In this respect, we confirm that 
the interaction between free excitations and disorder in topologically disordered systems (as long as ERMs describe 
them) has a peculiar mathematical structure that is different from disordered lattice systems (for lattice systems go 
vanishes identically) . 

To show this we shall compute the self-energ y pe rturbatively within two unrelated approaches: a) an improved form 
of the combinatorial formalism introduced in [la ], and b) a field-theoretic formulation. The field theory introduced 
here is quite different from standard formalisms in the theory of Random Matrices (see e.g. Refs. pliGjl]). It probably 
deserves an indepth study, which is left for future work. We remark that our combinatorial formalism is simpler than 
the field theory, and is probably the method of choice to carry out higher-order computations in the l/p expansion. 
However, it has the drawback that the asymptotic go(p 2 ) ~ P 4 appears at order l/p 2 from an exact cancellation of 
two contributions of order p 2 (at order l/p 3 we find an exact cancellation of ten contributions of order p 2 ). The 
field-theoretic framework clarifies that these cancellations are not accidental, and thus not restricted to low orders in 
the l/p expansion. 

The layout of the remaining part of this work is as follows: in sec. [IT] we discuss a particular phenomenon (phonons 
in topologically disordered systems) where a theory based on ERMs has been proposed in recent years. In sec. IIIII 
we anticipate our main result, namely the leading order of LmI7(p, z + i0 + ). In sec. IIVI we discuss in detail the 
combinatorial formalism up order l/p 2 . We describe the rules to group all the diagrams that arise at this order in 
a very small number of diagrams, according to their topological structure, and show that up to second order in the 
function go(p 2 ) the prefactor of the term oc p 2 cancels out. We also see that this cancellation appears in a given class 
of diagrams at l/p 3 . In order to shed a light over the mathematical origin of such cancellation, in sec. [V] we introduce 
a field-theoretical formulation that, despite producing a much larger number of diagrams, allows to give an argument 
explaining the origin of the cancellation at any perturbative order. 

II. A CASE STUDY FOR ERM: PHONONS IN TOPOLOGICALLY DISORDERED SYSTEMS 

Although ERMs have a wide range of application, in this paper we are mainly interested in the study of phonons 
in amorphous systems, such as glasses or supercooled liquids [24| . since the big amount of experimental evidences 
may provide fundamental insights about the correctness of the theory. Of particular interest is the case where the 
frequencies ui(jp) of the phonons with wave vector p lie in the GHz to the THz region (high-frequency sound). This is in 
fact the range explored by neutron and X-ray inelastic scattering experiments. These give the inelastic contribution 
to the dynamic structure factor, i.e. a Brillouin-likc peak with position to(p) and width r(p). Summarizing the 
experimental findings, for p < po (po is the first maximum of the static structure factor, typically a few nm _1 (25j ) 
one finds a linear dispersion relation oj(p) ~ cp, where the speed of sound c is quite close to that obtained by acoustic 
measurements. The dispersion relation typically saturates at p ~ po- Moreover, the p-dependence of the peak width 
is often described by r(p) oc p a . Interestingly enough, r(p) also saturates as the momentum becomes p ~ po- There 
has been a hot debate among different experimental groups about the value of the exponent a [27j , some claiming 
a ~ 2, and some a ~ 4. There is now some consensus that in the region where r is independent of temperature 
(i.e. w(p) > 1 THz) one has a = 4, while at lower frequencies (the GHz region), where r has a strong temperature 
dependence, the experimental value is a = 2 [28| . 

A simple model of the high-frequency sound is afforded by scalar harmonic vibrations around a topologically 
disordered structure made of N oscillation centers Xj, placed with uniform probability on a volume y[3l|. Particle 
displacements ifi have an elastic energy 



oo 




(10) 



n=0 




(11) 
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where the matrix M. has the form Eq. ([2]) and /(x) is the spring constant connecting particles separated by the 
vector x. We assume that /(x) is spherically symmetric, so that /(p) = g{p 2 ), where g is a smooth function. In the 
framework of the one-phonon approximation, the inelastic dynamic structure factor is related to the resolvent via 

„2 

S(p,u) = -— ImG(pV + iO + ). (12) 

UJTT 

As a consequence, the width of the Brillouin peak is related to the imaginary part of £ by 

lmZ(p,u>(p))=u>(p)r(p). (13) 

Then Eq. (|10[) implies that r(p) ~ p for very small p (for p ~ po the width saturates and a mixed, more com- 
plex scaling should be expected). In that regime the pho non-disorder interaction can be thought of as a scattering 
phenomenon of the Rayleigh type. 

Since ERMs describe the dynamics of vibrating particles within the context of the harmonic approximation, the 
theoretical predictions based on ERM theory must be compared with experiments in the region where r is independent 
of temperature; in fact the temperature dependence is an indication that the width of the peak is rather due to thermal 
processes, such as anharmonicities or relaxations, which require more refined theoretical approaches. 

We finally mention that vibrational frequencies u> are related to ERMs eigenvalues A (z = A + i0 + ), by the relation 
A = ui 2 , see Eq. (|T2|) . Hence, the width of spectral peaks in A-space and in w-space are related by Eq. (|13|) . 
Furthermore, Eqs. (j4|) and (fT0|) imply that the DOS in A space behaves for small A as g\(\) cx A' 13-2 '/ 2 , which 
translates to frequency space as a Debye spectrum g^{u>) cx uj ^ 1 (because of the Jacobian in the change of variable: 
dA = 2ujduj). At this point, the reader may object that lattice systems have a Debye spectrum even if go in Eq. (|10[) 
vanishes for them. In fact, their Debye spectrum is possible because Eq. ((4]) does not hold in the lattice case. 

III. THE MAIN RESULT 

The main result of this work is the following. Expanding the self-energy in powers of 1/p, i.e. 

S(jp,z) = S^(p,z) + S^(p,z) + ..., (14) 
where is of order l/p k , one has only one first-order contribution, 



I7«(p, z) = - f V(q, p)G (q, z)V(q, p), 

PJ (2tt) u 



(15) 



while at second order 

( P) z) = S f (P) z) + 4 2 ) (pj z) + 4 2 ) (pj z) > (16) 
where the three topologically different pieces are 



(17a) 



4 2 W) = - / ^^V(q,p)Go(q,z)IjM( qi z)G (c l ,z)V(q,p), 
PJ (2tt) u 

S b\p> z ) = -^ J -^^^p-^(p-q 5 P)G , o(q,^)^(q-k,q)G (k,z)^(p-q,p-q + k)x 

G (p-q + k,z)F(q-k,p), (17b) 
4 2) (p,z) = -i J ^^AAy( p _ q;P ) Go ( q;Z ) t /( q _ kj2q - p )Go(k,z)F(p-k,p). (17c) 



In Eqs. (IT5|) - (|T7|) we have used 

V(q,p)=p[f(q)-f(p-q)}, (18) 

which, as we will see below, plays the role of the interaction vertex. The bare propagator Go was defined in Eq. ^ . 
Note that V(q,p) ^ V(p, q). Other useful identities are 

V{q,p) = V{-q,-p), %p)=-y(p-q,p). (19) 



5 



Note that since V(q, 0) = 0, we have 

4 2) (0,z) = 4 2) (0,z) = 4 2) (0,z) = 0. (20) 

The high-density expansion for Laplacian ERM was introduced in [lH m. E q . was already reported there but, 
instead of Eq. (|16|) . one had 39 diagrams of order l/p 2 . Even if the final expressions were cumbersome, a numerical 
evaluation of the amplitude A in Eq. ([5]) was attempted for a simple choice of the function /. Presumably because of a 
numerical mistake, it was wrongly concluded that A ^= 0. Afterwards, it was announced (without supporting technical 
details) that the 39 diagrams previously found at order l/p 2 could be grouped as in Eq. (Q~6|) [18j |. Unfortunately, a 
numerical reevaluation of the amplitude A was not attempted from these simpler expressions. 

We remark as well that an independent computation of £ to order l/p 2 has appeared recently [3C|- We have 
checked that their results are consistent with ours, letting aside contact terms (in fact, these authors explictly state 
that some contact terms are lacking from their final expressions). Thus, their failure in identifying the g$ term in 
Eq. (|10p is not due to discrepancies in the final expressions. The underlying reason is rather more mundane, as we 
explain below. 

At first order the theory has the following behaviour. For small A, z = A + i0 + , we approximate the imaginary part 
of the propagator Go by 

ImG o (q,A + i0+) = -^^-^) (21) 

(assuming a linear dispersion relation e(p) ~ c 2 p 2 ). Then the only contribution to the imaginary part comes from 
q = VA/c. To evaluate the vertex V(q : p) at small q, and small p we just need to recall that /(p) = g(p 2 )- It is 
important to avoid any assumptions about the ratio p/q, which can be either very large or very small when both p 
and q are small (at the Brillouin peak p/q ~ 1, but in Ref. [30| it was unjustifiedly assumed that p<Cg). Then 

V(q,p) = g(q 2 ) - g(q 2 + p 2 -2p-q) 

« g(0) + g'(0)q 2 - g(0) - g'(0)[q 2 + P 2 - 2p ■ q] 

= -g'(0)[p 2 -2p-q]. (22) 

If we now square the vertex function and perform the angular integral, we obtain (Sd is the surface of the sphere in 
D dimensions) 

[sW^lV+i? 2 /)- (23) 



D 

The integral over q is now straightforward, thanks to Dirac's S function in Eq. (fT5|) . We get 



Imi; (1) (p, A + i0+) 



(24) 



Hence, already at order l/p, go(p 2 ) in Eq. (|10l) is of order p A . Had we neglected the p 4 term in confront of the q 2 p 2 
term (as done in Ref. [30j). we would have failed in identyfing the go term. The physical reason for which the presence 
of such a term is mandatory (namely the existence of a Debye spectrum), was discussed in the concluding paragraph 
of Sect.HH 

Let us now check that the amplitude A in Eq. (|HJ| vanishes. We merely need to compute the imaginary part of the 
self-energy at its lowest order in A, namely \~ ~ . A general term of the diagrammatic expansion involves the factor 



d u qG (q,z) = jdqdf2 D q n - L G (q,z), (25) 

where we have expressed the measure in terms of polar coordinates in D dimensions. Every bare propagator is 
associated to one or more vertices that are smooth functions of the involved momenta. In fact we can expand the 
product of such vertices in a Taylor series. Now, the point is that if we want the lowest order in z, we have to exclude 
all the terms that are proportional to q and we have to take only the zeroth order term of the Taylor expansion. We 
can obtain this term simply making the following substitution 



Im J d C qG o (q,A + i0 + ) = -1\<P-*V* J dqdQ L 



■S(q) (26) 
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whore Im stands for the imaginary part proportional to X [D 2,//2 . Then 

Im4 2) (p,A + iO+) = -|a^~ 2 )/ 2 / dq J df2 D V 2 (q,p) oc X^^p 2 



and 



lmS {2) (p, A + i0+) = -ImE%> (p, A + 10+) . 



(2) 

It follows that the amplitude A vanishes, because the U c contribution is already of order p 

lmS {2) (p, A + i0+) oc -2A^- 2 )/ V(p, p)^ 1 ' (p, z) oc \( D ~ 2 ^ 2 p i , 



(27) 



(28) 



(29) 



In the following, we will show explicitily that the cancellation of the \( D ~ 2 ^ 2 p 2 term arises even for a given topo- 
logical class (quite large) of diagrams at 1/p 3 order, and we will provide an argument that predicts such cancellation 
at any pertubative order. 

IV. THE COMBINATORIAL COMPUTATION 

The first approach to the computation of the resolvent is based on the expansion of Eq. ([3]) as a power series, 



N / N 

G (p> *) = E psr n E eip ' (Xj " Xi) w 

R=l \ i,j=l 



(30) 



Although the final results will only depend on p, in order to develop the formalism it is convenient to reintroduce the 
dependence on p. 



A. Organizing the calculation. The bare propagator. 

1. Momentum shift: choosing wisely the integration order 
The i?-th term of the expansion Eq. (l30l) is 



1 N 

E TV E e ip '^«- x *o) [M i0 , il M il , ia ...Mi R _ 1 ,i R \, (31) 

i ,i 1 ,...i R i,j=l 

where the average over the vibrational centers take the form of a multi-dimensional integral with measure 

i r N 

J i=l 

As for all such integrals, although the final result is independent of the order in which the individual integrals are 
performed, the difficulties encountered in a real computation are dramatically smaller if one finds a wise ordering for 
iterated integrations. 

Now consider the expression 



(32) 



which arises as a factor when introducing the explicit form (Eq. [2]) of Ai into Eq. (|3"Tj) . When dealing with a diagonal 
term, we shall integrate over the position of the medium particle, X&, ; when dealing with an off-diagonal term, we shall 
integrate over Xj, . For a diagonal term, the integral over the position of the medium particle is easy, if the particle 
index ki does not appear elsewhere in the chain (even if the index i\ is sure to appear at least once more along the 
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chain). For the non-diagonal term, the integral over x^ ; is very simple if it does not appear later in the chain (even if 
appears twice or more times in the chain, to the right). The two integrals yield 



V 



[f(0)5^ l+1 -[l-6 liai+1 ]f(p)] C -^ + 



Since a term of order R has R such factors, the number of values the index fc; (or i{) can take without violating the 
non-repetition condition is between N and N — R. But both N/V and (N — R)/V tend to p in the thermodynamic 
limit, hence momentum can shift through non-index-repcating elements from left-to-right: 



-ip-Xi 



P[f(0) - /(p)]e 



(33) 



Similarly, momentum can shift through non-repeating elements from right-to- left. In that case, one would integrate 
over Xfc ; (diagonal term) or over Xj, (non-diagonal): 



e ip - x ^+i 



«p[/(0)-/(p)] 



(34) 



Note that a given matrix-element might be considered as non-repeating for momentum shift from right-to- left, but it 
could be not suitable for the left-to-right momentum shift. 

At this point, the computation of the leading order is straightforward. If there are no obstacles for momentum 
shift, we just push to the left exponential c~ lpXi o to the right until it cancels out with e lp Xi «, leaving us with (since 
there are precisely R matrix elements) 



p R [f(0)-f(p)] R . 
Then the high-density limit of the sum in Eq. ([30]) is 

1 



G (p,z) 



z - e(p) ' 

which is then the bare propagator of the theory, as anounced in Eq. (jB]). 



(35) 



(36) 



2. Repeated indices 

Now consider a situation where we can shift the external momentum p from left-to-right until a particular particle- 
index (say %i = 1 or ki = 1) is repeated in the chain somewhere to the right, so that we must stop. At this point, we 
shift the external momentum from right-to- left, until a particle-label repetition i r +i = 2 or k r = 2 stop us. We depict 
this situation as 



1 [stuff] 2 . 



(37) 



At this point we will have 



p s [f(0) - /(p)] s e-P' x n 1 [stuff] 2e ip - x '" +1 /[/(O) - /(p)] J 



Since the very same scheme of particle-label repetitions 1 [stuff] 2 can be found for all values L, S — 0, 1, 2, . . ., we can 
sum all those terms to find a contribution 

G (p, z)c- ip X! ' 1 [stuff] 2c ip - x *-+i G (p, z). 

We interpret the two factors Go(p, z) as the external legs for a Dyson resummation of the self-energy. 

Clearly particle-label repetitions are going to be very important in what follows, so some terminology will be useful. 
A generic factor 



J/(Xi 



J 



(38) 
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will be called an L-stop if particles ki or i\ are repeated somewhere to the right (so that momentum cannot be shifted 
from left-to-right trough index i\). Similarly, we shall call it an R-stop if fc; or are repeated somewhere to the 
left. We note that a matrix element can be both a L-stop and a R-stop (if ki is repeated both to the right and to the 
left, or if ii is repeated to the right while is repeated to the left). 

To make momentum flow through an L-stop or R-stops we resort to the so called fake-particle trick. Consider a 
particle label, say 1, that appears twice (for instance, in an L-stop and in an R-stop to its right). Before carrying out 
the average over {x^}, we multiply the term by 1 written as 

1 = J A D yi 5( Xl - yj) = j d'Vid^qe^"!). (39) 

Then we can pretend that particle 1 takes two identities, 1 and 1, so that there is no repetition. The price we pay for 
this simplification is that: 

• we have an extra integration over the internal momentum q, 

• we have to deal with an extra factor e lqxi at the L-stop, and an extra e _lqy i at the R-stop, and 

• the fake particle yj does not bring a combinatorial TV factor, or an l/V from the normalization of the y 1 integral, 
so that there is a lacking factor of p (this we can ignore if we add compensating 1/p to the final expression). 

However, the modified momentum-shift formulae are simple enough to justify these inconveniences. Integrating over 
xi we obtain 

e-ipx^iqx, [^/(xx-ij,) - (l-£ Ml+ J£ Ml /(xi-^ +1 )] — ► V(q,p)e- i ( p - q K +1 . (40a) 
Similarly, integrating over yj at the R-stop, we have 

£v,v +1 /(yi ~ x v) - C 1 - S i,ir)k,i r+ J( x ir -yi) 

As a warning on momentum-shift, note that one may shift momentum from left to right as long as there is nothing 
to the left still needing integration (and similarly for right-to- left shifts). Momentum shift can be visualized as a zip 
with two heads: one pulls both heads until they meet (and then there are no integrals left to be done). 

3. The reduction formula 

Imagine we face the situation 

...l[stuff]l..., 

i.e. the leftmost stop is an L-stop at index position I, the rightmost stop is an R-stop at index position r + 1, and 
the particle that prevents the two momentum shifts is the same at both ends, say ii = i r +\ = 1 or any other possible 
combination (fc; = i r +\ = 1, %i = k r = 1, or fc; = k r ). If the particle label 1 does not appear again inside the brackets, 
a nice reduction formula follows: 

c"' px,i [^,i, +1 /(xi -x i; ) - [l-5 lji(+1 ]5 i|) i/(xi -x i(+1 )] x [stuff]x 

[<5 ir , ir+I /( Xl - Xi r ) - [1 - ^ijS^i/K - xi)] e ipx v+i _> 

~p /(0j^(q.P)^ iqXl ' +1 Nuff]c +i ^. (41) 

This can be proved by averaging over Xi. To adjust the power of p, just recall that there were order TV choices for 
(say) the index coincidence fc; = i r+1 = 1. Note that the proof of Eq. (|4"T|) involves doing four different integrals. 
Let us see how the fake-particle formulae (Eqs. (j40])) yield the same result effortlessly. The introduction of the fake 
particle transforms the left-hand side of Eq. (|4"Tj) to 

c-ipx^iqx: [ SiiA+ J( Xl - x. i; ) - (1 - <5 Mi+1 )W(xi - x [stuff] x 

Si r ,ir+J(yi - X ir) ~ i 1 ~ S 1 ,i r) S l ,l r+ 1 I ( X V ~ Yl) 



c ipx Ir+1 c -iq yi _^ e i(p-q)x ir p j ( 4Qb ) 



e ipx Ir+le -iq yi ( 42 ) 



9 



We now merely shift momentum to the right using Eq. (|40a[) and to the left using Eq. (|40b[) to obtain 

^l/(q,p) 2 e" i(p " q)Xl '+i x [stuff] x c +Hp-q)x^ 

Now a change of integration variable q — > p q and the second of identities (fT!?]) yield Eq. (|4T|) . 

Both EM and E^ follow directly from Eq. $2]). Also, more general expresions can be found easily from it, as we 
shall see below. 




B. Order 1/p 

If no further particle-label repetition arise, the momentum e _lq Xil + 1 in Eq. (|41[) can be shifted to the right until it 
is killed by the second exponential e lq . We have then a set of contributions of the form 

v 1M/(Q) - Ttt x i [g*^ p) 1AM -fmr x W(°) - /(p))1- _ (43) 

z a+i oJ(2tt) d v ' z b+1 z c+1 v ' 

a+b+c+2=R r J \ ) 

composed of the product of three harmonic series that are easily seen to add-up to 

G (p, *)*- J- T^p^fa' ?) G ^ z ) x G °(P' z )- ( 44 ) 
If we interpret the two factors Go(p, z) as external legs of a Dyson resummation, we get 



E^(p,z) = - /AAy 2 (q lP )Go(q, 



z), (45) 



which is the first order result anticipated in sec. 



C. Order 1/p 2 



The cases with two pairs of repeated indices, or one index occurring three times contribute to the second-order 
corrections. The contributions separate naturally in three kinds, according to the arrangement of the repeated indices. 



(2) 

1. The nested case: S A 

Take now the scheme of particle repetitions giving rise to Eq. (|43j) . and place it in between an external pair of 
particle repetitions: 

...2...1...1...2.... 

Assume that the index 2 happens twice and only twice in the chain. We are thus entitled to use the reduction formula, 
Eq. (|4"Tj) , for particle 2. The inner momentum q, can then be shifted (from cither side) until it hits particle 1, where 
it produces a contribution such as Eq. (|4"5)l . The only difference is in that the role previously played by the external 
momentum p is now played by the internal momentum q. We get, without need for further computation, 

Ef{ V ,z) = - f^V(q,p)Go(q,z)E^\q,z)G (q,z)V(q,p). (46) 
9 J (271T 



2. The interleaved case: E B 

The E^ piece in Eq. (|16p arises from the pattern 

...I...2...I...2.... 
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A moment's thought indicates that the leftmost 1 must belong to an L-stop, while the rightmost 2 must be an R-stop. 
Furthermore, the internal 2 should be an L-stop (otherwise, one would use a fake particle to shift momentum from 
left-to right over it trivially). For the same reason, the internal 1 should belong to an R-stop. 

Our previous succes with the reduction formula, Eq. (|41[) suggests that we try to deal with all such terms at once, 
by performing the integral 



d^ Xl d D x 2 _. 

e 1 



V 2 



[<^a+i/( x i - x iJ - U - <*i,i I+ i)£i,i,/(xi - [. . .] x 

[£v,i r+1 /(x 2 - x ir ) - (1 - 5 2 , tr+1 )5 2 . ir f(x 2 - x Jr+1 )] [. . .] x 
[^i s ,i, + i/(xi -Xi„) - (1 - Si,i s )Si,i a+1 f(xi a - xi)] [. . .] x 

[^ iWl /(x 2 - Xi J - (l-<y2,i,)*2,i, +1 /(xi, -x 2 )] c ipXi *+i. (47) 

Here the several [. . .] stand for arbitrary numbers of matrix elements without momentum stops arises. Note that not 
all the ...1...2...1...2... terms have the form Eq. (|4"7]) : the central 2, 1 particles may also collapse onto a single 
matrix clement (necessarily non-diagonal) which is both an R-stop and an L-stop. Such terms will be considered in 
sec. IIVC 31 

We now introduce two fake particles, 1 and 2, to transform the above integral into 



d^qd^k /"d^xid^yjd^xsd^ya 



(27T) 



2D 



K,,: 1+1 e iqxi /(xi - x 4; ) - (1 - <5 lli!+1 )<5i )i( e^/(xi - [...] x 



e 4fcX2 /(x 2 - x 4 J - (1 - 5 2 . tr+1 )6 2 , lr e^f(x 2 - x* r+1 )] [. 
^-^/(yj - xj - (1 - Jj^)*! e-^/(x is - yi )l [. 



x 
1 x 



f(y 2 - xij - (1 - <y 2 , i ,)^, i . +1 e- <kys /(xi. - y 2 ) 



One then shifts momentum from left-to-right up to i s and from right-to-left again up to i s , to find 
1 /"d D qd D k 



V(q, p)G (p - q, z)V(k, p - q)G (p - q - k, z)f (q, p - k)G (p - k, z)F(k, p) 



(27r) 2i5 

Eq. (|17b|) follows after changing integration variables according to 

p — q — > q, and q k 



(48) 



(49) 



(50) 



3. T/ie collapse of a L-stop and a R-stop: S c 

As we have remarked, it can happen that the L-stop and R-stop of Eq. (|47p actually belong to the same matrix 
element, necessarily non-diagonal. However, any non-diagonal term should be paired with a diagonal one. As we 
mentioned in sec. IIV A II a diagonal term can be both a L-stop and a R-stop if the medium particle is repeated both 
to the left and to the right. Hence we will be considering here this kind of terms (0:of diagonal matrix element, D: 
diagonal matrix element): 

1... 0(21)... 2 + 1...L>(1)...1 (51) 

The the terms with an off-diagonal index appearing three times (1 . . . 0(1?) ... 1) do not belong to S c , and are 
considered in sec. IIV C 41 

Let us start with the case 1 . . . D(l) ... 1: 

J — — e~ ipx '' [<5. Ml+1 /(xi -x^) - (1 - <5i iil+1 )5i >i! /(xi - Xi I+1 )] [. . .] x 

Si r ,i r+1 f(xi - x ir ) [...] x 

[<$w, +1 /(xi - x ls ) - (1 - <5 M J<5 Ms+1 /(x ls - xi)] e ipXi »+i. (52) 

We now introduce two extra fake particles to substitute particle 1, namely 1 and 1 via the identity 

1 = Jd D yi |d D z i 5(x 1 -y i )5(y i -z i ) = ^p Jd D qd D kd D yi d D z i ^-y^k^-z^ (53) 
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to find 



Ji0^ J * Zl e " iPXi ' fo-^K* ^ <^ (+1 )<W(*i - x, +1 )] e-» [. . .] 



i(k— q)xi 



-ikz.~ 



(54) 



^ 8 ,i s+ i/(zi - x ls ) - (1 - ^ jis+i )5i jis+i /(x is - zj) 
Finally we shift momentum from lcft-to right up to i r , from right-to-left up to i r +i and integrate over Xi. We obtain 

1 fd D qd D k^ 



l...D(l)...l = -, 



P 2 J (27T) 



2/) 



T(q, p)G (p - q, z)[pf(k - q)]G (p - k)F(k, p) 



Consider now 1 . . . 0(21) ... 2: 
d D xid D x 2 



Shifting momentum left-to-right up to i r and from right-to- left up to i r +i we are left with 

- 1 fd D x 2 d D yi /(x 2 - y^e-CP-d-^Cx.-xO = _^ (p q k) 



so that adjusting the power of p, we get 

1 fd D qd D k 



1... 0(21)... 2 



(27T) 



2L> 



F(q, p)G (p - q, z)[-p/(p q k)]G (p - k, z)V(k, p). 



Adding together the two pieces, Eqs. (|59p and ([55]) we finally find 
l(2 )_ 1 fd D qd D k 



s, 



I ( 27r )2D P)°o(P ~ q, «)V(k - q, p 2q)G (p - k, z)V(k, p), 



which after the change of variables 



and use of identities Eq. (fl"9"|) yield Eq. (|17c|l . 



p - q, 



p k, 



(55) 



y2 - " ipXli [<5 liA+1 /( Xl - x i; ) - (1 - ^i, i(+1 )5i, i( /(x 1 - Xi l+ J] [. . .] x 

[-<5 v , 2 R +1 ,i/(x 2 - Xl ) [. . .] [<5 is , is+1 /(x 2 - x is ) - (1 - 5 2 , is )5 2 , is+1 /(x is - x 2 )] c ipx< -+i. (56) 
Introducing two fake particles, 1 and 2 we can rewrite it as 

J^^fDV* * 2 6 " iPXi ' elqXl [^A+^( Xl - Xii) - (1 - <W)<W(*i - *i, + J] [•••]>< 

[-e" iq ^] 5v >2 5 ir+1 ,i/(x 2 - yi )e ik ^ [. . .] k,i. +1 /(x 2 - - (1 - &s,i>3,i. +1 /(*i. - y 2 )j e iky *e** • . (57) 



(58) 



(59) 



(60) 



(61) 



4- The Dyson resummtion to order l/p 2 

Recalling Eq. (fT4"]) , we notice that we have still not identified the pattern of particle-label repetitions that gives rise 
to the second-order terms appearing in the Dyson resummation the first-order self-energy, 



£«(p,z)Go( P ,z)£ (1 W), 
(we have not written the irrelevant external legs). The natural candidate is 

1 ... 1 ... 2 ... 2, 



(62) 



(63) 



where the sequence is L-stop, R-stop, L-stop, R-stop. This expectation is correct, but it will turn out that the 
constraint imposed by the matrix-product structure needs extra terms to build Eq. (|62p . These missing terms will be 
provided by the pattern 1 . . . 0(1?) ... 1. 
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Let us first compute blindly the term 1 ... 1 ... 2 ... 2, incurring in a quite instructive mistake. We introduce only 
one fake particle, 1: 

e -ipx ile i qxi [5. ii|+i/ ( Xl _ x . ti ) _ (i _ 5 1A+1 )6 1>it f( Xl x il+1 )] [. . .] x 

e" iqyi [Si^JiYi - *ir) - (1 - h,i r )k,i r+ J(xir - Yl)} [■■■}>< 
[&„i„+i/( X 2 - XjJ - (1 - S2 t i s+1 )S 2 ,i s f(X2 - x is + 1 )] [...] x 

[8i a , it+ J(x2 - x ir ) - (1 - fc.O&.u+i/fo - x 2 )] c ipx ^+i. (64) 

We shift momentum from left to right up to i r as usual. At this point, we still need to push the momentum to the 
right (this is unusual). We need to perform two integrals, 



d°x, e 



-i(p-q)x; -iqyj 



Hence the integrations up to this point yield 



-(i - k, ir )k,i r+ J(xir yd] = -Hp q)e ipx -+i 



p J (2tt) u 



(65) 
(66) 

(67) 



It seems to be an easy matter to complete the computation: one pushes momentum p to the right up to i s , seemingly 
yielding a bare propagator Go(p, z), and we would be left with 2 ... 2 (a standard diagram for the self-energy at order 
1/p). However, after some reflection it is clear that an R-stop and an L-stop such as . . . 1 . . . 2 . . ., where both particle 
1 and particle 2 appear in off- diagonal matrix elements must be separated by at least one off-diagonal matrix clement. 
Hence if there are S matrix elements between the R-stop and the L-stop, when shifting momentum p we will encounter 
a factor 



p s [f(0)-f(p)} s -[pf(0)} 



which adding the geometric series means 



G {p,z) 



z - Pf(0) 



Hence the correct result is 

1 ... 1 ... 2 ... 2 = sW(p, z)G (p, z)S^(p, z) 
1 1 /YPqd^k 



, y(q,p)Go(p-q^)[p/(p-q)]b/(p-k)]Go(p-k,^y(k,p). (68) 

We will now show that the second term in Eq. (|68[) is cancelled by the contribution from 

1...0(1?)...1. 

In this pattern, the leftmost 1 belongs to an L-stop and the rightmost one to an R-stop. The first observation is 
that the central 1 in the 0(1?) must nessarily appear in an R-stop (because we never find the same particle in any 
matrix element /(xj — Xj), and due to the constraint imposed by the matrix product). The second observation is 
that there must be, at least, one off-diagonal matrix element between the two R-stops sharing the common particle 
1. Introducing fake particles 1 and 1, we are left with 



d^qd^k fd^Xjd^Vjd^zj 
] V 



c- ipx " [<fi„i 1+1 /(xi -*,) - (1 - 5 Mi+1 )<W(*i - e iqXl [. ■ •] 



-(i-^ r+1 )^ >ir+1 /(xv-yi) 



3 i(k-q)yj 



ikz=; 



s+ i/( z i - XiJ - (l - h,i 3 ) s i,i s+ J( x is - z i) 



lpXj 



(69) 
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All that remains is a simple momentum shift, keeping in mind that when going over the factor [...]' it will give 

G (p-k,z) L^=- p /(p-k) (?o(p ~ k ' z) . (70) 

z-pf(0) m z-pf(0) 

Thus one finally finds 

1 1 /"d^nd^k 

1 . . . 0(1?) . . . 1 = j—- / q U(q, p)G (p - q, z)[pf(p - q)][[p/(p - k)]G (p - k, z)V(k, p). (71) 

z — pf(0)P J 



D. Higher orders 



The argument of sec. IIV C II is fully general. Consider the contribution of order l/p n to the propagator, rather 
than the self-energy (i.e. let us include both the connected and disconnected pieces). We can write this as 
Go(p, z)H^ n \p, z)Go(p, z). Let us emphasize that H^(p,z) refers to the full contribution to the propagator at 
order 1/p", not to a particular topological subset (such as the cactus (l7|). 

We may enclose the scheme of particle label repetitions that generates H W (p, z) within an L-stop and an R-stop 
with equal particle labels that do not appear again along the chain. Under such circumstances, we are entitled to use 
the reduction formula, Eq. (|41[) . which yields 

4" +1) (p,z) = - / 7 ^^(q,p)Go(q,z)iJ ( " ) (q,z)G (q,z)y(q,p). (72) 
pj (2ir) u 

Clearly this is not the full self-energy at order l/p n+1 , but it is a genuine part of it that automatically verifies 

4" +1) (p,z) = 0. (73) 
In particular if n = 1 this gives the 1/p 2 term S\(p, z) discussed above. It is interesting to note that 

Im4" +1) (p, A + i0+) ~ p 2 X^ D ~ 2 y\ (74) 

since for q ~ 1, and z = A + i0 + , for small A it is expected that (Debye spectrum, see Sect. [TTJ) 

MW(q,z)«A( fl - 2 )/ 2 . (75) 

Thus, the vanishing of the amplitude A in Eq. ([5J implies that non trivial cancellations occur at all orders in 
perturbation theory. Since we are presenting arguments for such cancellation, we agree with Ref. [3(J in that the 

ter ms alone do not reproduce the correct analytic structure of the theory. 

1. Towards the self-energy at third order 

Using the combinatorial rules described above, it is possible to push the perturbativc computation to order 1/p 3 
order, which has never been attempted before. Here we will limit ourselves to the terms without collapse of an i?-stop 
with an L-stop (i.e. we will retain only the terms with 6 vertex functions). The reason is that the combinatorial com- 
putation suggests very simple Feynman rules that can be used to obtain the diagrams, without lengthy computations. 
The purpose is to check that, at least within this subclass of diagrams, the cancellation of the prefactor of the p 2 ui D ~ 2 
term still occurs. 

Let us describe the Feynman rules. Take for instance a term such as 

L1...L3...L2...R3...R1...R2. 

The rules are as follows 

1. Draw an horizontal full line and mark on it the stops (preserving the ordering). 

2. Join the corresponding L- and R-stops with a dashed line. 

3. The diagram has an incoming momentum p, from left to right. 
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4. Attach a momentum to every line (full or dashed), applying momentum conservation at each stop. 

5. Associate a bare propagator, Go, to each full line. 

6. Associate a vertex function to every stop, such that its first argument is always the momentum runing over the 
dotted line. 

7. For an L-stop, the second argument of the vertex will be the momentum running over the full line to its left. 

8. For an R-stop, the second argument of the vertex will be the momentum running over the full line to its right. 

9. Multiply by 1/p 3 and integrate over the internal momenta. 

Applying these rules to the patterns without stop collapse, we obtain the following contributions. 
a. Terms LI . . . L2 . . . L3 . . . R3 . . . R2 . . . Rl 



\ \ | i rr\ D ac\ D \<:c\ D \ 

= ' ' ' =^3 J ,y V (p-q,p)g (q,z)V(q-k.q)go(k,z)y(k-l,k)x 

G (l, z)V(k - 1, k)G (k, z)V(q - k, q)G (q, z)V(p - q, p). (76) 



Now we wish to compute (Sd ■ surface of the unit-sphere in D dimensions) the limit 



(2n) D lmh(p,z) 

Jl = -^s c 7}^^U^W (77) 



and in general, defined from Ik(p, z) as the same limiting procedure. 
The rules to obtain the limit painlessly are simple: 



1. Locate a propagator, G{q) whose running momentum is never a second argument of a vertex function V(-,q). 

■kSd , 
' {2tt)D' 



2. Substitute that propagator by — Z P D z^ D 2 ^ 2 5(q), and perform the q integral. 



3. Apply the simplification 



G (q, ^(q, q) = ^ q) = 1 + 0{z) . (78) 



For 1\ only I = gives a contribution to Ji, hence 



Jl = 7 I d ( 2 ^)^ k F(P - q ' p)Go(q ' z)F(q - k ' q)F(q - k ' q)Go(q ' z)F(p - q ' P) ' (79) 



b. Terms LI . . . L2 . . . L3 . . . R2 . . . R3 . . . Rl 



/ / / x v 



; / i s s v 

h = ' 1 = ^ J - ( q 27r)3J V(p - q, p)G (q, z)V(q - k, q)G (k, z)V(k - l,k)x 

G (l, z)V(q - k, q - k + l)G (q - k + 1, z)V(k - 1, q)G (q, z)V(p - q, p). (80) 
For I2 one easily realizes that only I — contributes to J2: 

J ' 2 = ^J (2^p F(P ~ q ' P)G ° (q ' Z) ^ (q " k ' q) ^ (k ' q)Go(q ' Z)V(P ~ q ' P) (81) 
Since V^(q — k, q) = — V(k, q), one has J2 = — Ji. 
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c. Terms LI . . . L2 . . . R2 . . . L3 . . . R3 . . . Rl 



f I 1 /"d^qd^kd^l 



/d ad kd 1 
(27r)3J V^P - q, p)G (q, z)V(q - k, q)G (k, z)V(q - k,q)x 

G (q, z)V(q - 1, q)G (l, z)V(q - 1, q)G (q, z)F(p - q, p). (82) 



Both k = and 1 = contribute to J3 (for k = we changed the dummy variable 1 to k) 



2 /cPad^k 

■h = -j J ^ )2D V(p - q, p)V(q - k, q)G (k, z)V(q - k, q)G (q, z)V(p - q, p). (83) 



d. Terms LI . . . L3 . . . R3 . . . L2 . . . Rl . . . R2 



1 — ^ / 1 N 1 M ad kd 1 

h=\ 1— i 1 ! * =-3 y (2?r)3I3 F(p-q, P )Go(q,z)F(q-k,q)G (k^)F(q-k,q)x 

G (q, z)V(q - 1, q)G (l, z)V(p - q, p - q + l)G (p q + 1, z)^(q - I, p). (84) 



For J4 both k = and 1 = arc relevant 
W~ ad 



1 fd qd k 

1/4 = ^3 y (2^pJ ^ P ~~ q ' p ) y ( q ~ k ' ^^(k, z)V(p q, p q + k)G (p - q + k, z)V(q - k, p) 



1 rd D a d D k 

-3 y (2 ^ )2J V(p - q, p)G (q, z)V(q - k, q)G (k, z)V(q - k, q)V(q, p). (85) 



e. Terms LI . . . L2 . . . L3 . . . R3 . . . Rl . . . R2 
" ~ *C 

/ s 
/ / 
/ t 



\ X 
\ \ 
\ \ 



h=[ [ CZl ! ! =j 3 J - ™ - F(p-q,p)Go(q,z)F(q-k,q)G (k^)F(k-l,k)x 

G (l, z)V(k - 1, k)G (k, z)V(p - q, p - q + k)G (p - q + k, z)V(q - k, p). (86) 



The only relevant contributions is now 1 = 



fd a d k 

' h = 7 j (2tt) 2 ^ V{P ~ q ' P)G ° (q ' Z)V{<1 ~ k ' q) ^ (P q ' P q + k)G ° (P q + k ' Z)V{<1 KP) - (87) 
/. Terms LI . . . L2 . . . Rl . . . L3 . . . R3 . . . R2 



1 I \ /-* x \ 1 /-d^qd^kd^l 



J 6 =l 1 ! 1 i ! =^3 y (2^)3£> ^(p-q,P)go(q,z)V r (q-k,q)Go(k,z)y(p-q,p-q+k)x 

G (p - q + k, z)V(p - q + k-l,p-q + k)G (l, z)V(p - q+k-l,p-q + k)G (p - q + k, z)V(q - k, p). 



Jq stems both from k = and from 1 = 0. For the k = contribution, wc make the change of variable q — > p — q 
to identify the cancellation with J3: 



1 fd n a d^k 

J 6 = -j y { ^ )2D V(q, p)V(q - k, q)G (k, z)V(q - k, q)G (q, z)V(p - q, p) 
1 r d D qd 

7 J ~W) 



1 r d^od^k 

J (2*)™ ViP ~ q ' P)G ° (q ' Z)V(q ~ k ' q)G ° (k ' Z)y(P q ' P q + k)t/(q ~ k ' P) " (89) 



g. Terms LI . . . L2 . . . LZ . . . Rl . . . RZ . . . R2 



7 7 =, 1 1 ! 5 1 = ± / (2 ° )3J y(p-q,p)Go(q,g)V(q-k,q)G (k,z)y(k-l,k)x 

G (l, z)V(p - q, p q + l)G (p q + 1, z)y(k 1, p q + k)G (p q + k, z)V{q - k, p). 



Only 1 = contributes to J7: 

iD„ A D 



1 Ai qcPk 

■7? = -j y ^ )2g y(p - q, p)G (q, «)^(q - k, q)V(k, p q + k)G (p - q + k, z)y(q - k, p) 



h. Terms LI . . . LZ . . . L2 . . . Rl . . . RZ . . . R2 



/ / / 



XXX 



/ / / v v 1 /"d^nd^kd^l 

J 8 =l 1 1 ! ! 1 =-3 y y(p-q,p)G (q,z)y(q-k,q)G (k,z)^(k-l,k)x 

G (l, z)y(p q, p q + l)G (p q + 1, z)V(q k, p + 1 k)G (p + 1 - k, z)V(k - 1, p). 



Again, only 1 = matters: 

^D„ ad 



1 I'd qd k 

J 8 = -3 y ^ 2£ , V^(p - q, p)G (q, z)V(q - k, q)V(q - k,p - k)G (p - k, z)V(k, p) 



i. Terms LI . . . LZ . . . L2 . . . RZ . . . Rl . . . R2 



1 / / v n v v 1 /YI £l ad- D kd' D l 

/ 9 =l 1 1 ^ ! ! =-r y (27r) 3D y(P-q,P)Go(q^)V(q-k,q)G (k,2)7(k-l,k)x 

G (l, z)V(q - k, q k + l)G (q k + 1, z)V(p - q, p - k + l)G (p - k + 1, z)V(k - 1, p). 



And, once again, only 1 = contributes 

iD„ A D 



1 fd^nd^k 

Jo = -r y (2 ^ )2D ^(P - q, p)G (q, z)V(q - k, q)y(p q, p k)G (p - k, z)V(k, p) 



j. Terms LI . . . L3 . . . Rl . . . L2 . . . RZ . . . R2 



7 / x / \ x 2 (d qd kd 1 

ho=l 1 ^ 1 1 = J — ( 2 7r) 3D — F (P _t l'P) G o(q,2)V"(q-k,q)Go(k,z)^(p-q,p-q 

Go(p-q + k,z)y(p-q + k-l,p-q+k)Go(l,z)V(q-k,q-k + l)G (q-k + l,z)y(p-q + k-l,p) 



Here we have a contribution from k = as well as from 1 = 0: 



2 fd qd k 

Jl0 = J {2n) 2D V( - P ~ q ' P ^ P q k ' P ^o^' Z ) V (^ q + k ) G o(q + k, z)y(p - q k, p) 
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2. Resummation of the imaginary parts 

The resummation of the imaginary parts of the previous diagrams is simple. Using the properties of the vertex 
V(p, q) and changing carefully the integration variables when necessary we can show that 

Ji + J 2 = (98) 

J3 + J4 + J6 + Jw = (99) 

J 5 + h = (100) 

Js + J9 = (101) 



so that 



in 



£ Ji = 0, (102) 

i=l 

and the total contribution to the imaginary part proportional to z^ -2 )/ 2 vanishes. 



V. A FIELD THEORY APPROACH 



In this section we will introduce a field-theoretical representation for the resolvent G(p, z). Within this formalism 
one is able to obtain the perturbativc computation for the self-energy in a more strightforward way than with previous 
formulations fl5T ]. Interestingly enough, due to the ultraviolet behaviour of the bare propagator of the field involved, 
such perturbative expansion yields some divergent terms that can be summed up to zero. The starting point is the 
following representation for the resolvent: 



1 _ / [Xli d(j)i) fafa exp {-i J2hn </>l \( z ~ Efe /( x i - x fc)) <W + /(*( - Xm)] 0m} 

G(p,z) = -^X) e<P ' (x * 



Introducing the fields 



one has 



/ (ilf ex P {-5 Sim K Z - Efc /( x * - x fe)) + ./"( x i - x m)] </>m } 



^ (X)= {arbitrary x + x, 

p( x ) = -X^( x_xfc )' 



G(p,z) = ^ | d^xd^y e ip-(x-y) PWp(y) | |j]^( X?; )0(x)0(y) exp{5 p [0]} 



(103) 



(104) 
(105) 



(106) 



where we have introduced the action and the partition function at a fixed realization of the disorder, given respectively 
by 



s P W\ 



f /d c xd^v0(x) 



zp(x)5(x - y) - pp(x)S(x - y) j& D o- /(x - er)p(er) + pp(x)/(x - y)p(y) 



(107) 
(108) 



Now we note that the action Eq. f|10T[) depends on the field <j> only through the values that it assumes on the 
random positions {xj}. In fact, in the action, the field (f> is always multiplied by the random field p, which selects the 
random points of the lattice {x^}. So, we can substitute the discretized functional measure with the continuous one: 
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this is a crucial step. The continuous version of the fuctional integral is invariant under the following transformation 
of the field cf>, which we shall call gauge transformation: 



6'(x) = <f>(x) + h(x), with h(xi) = 0, i = l...N. 



(109) 



This is a local transformation, but we can see that its global version is trivial because the condition Eq. (|109[) implies 
that a global transformation can be possible only for h = 0. This local symmetry is not present in other field-theoretic 
formulations (lEj . 

We now look at the resolvent: it can be written in the form 



G(p,z) = fd D xd D ye^^ ^p(x)p(y) (0(x)0(y)), 



(110) 



where (•) stands for the average over the action S p [(f>]. We immediately see that p(x)p(y)((j)(x)(j)(y)') is gauge invariant. 
With the change of variables 



p(x) = 1 + 5p(x) 



the resolvent can be written 



G(p,z) = Jd D xd D ye^^ ^ (0(x)0(y)) + 2^Sp(x) (0(x)0(y)> + ^Sp(x)6p(y) (0(x)0(y)) 



(111) 



(112) 



with the action 



S P [<f>] 



d^xd^y 0(x) [(z - pf(0))6(x - y) + pf(x - y)j <j>(y)+ 
P - J d D xd D yda (/>(x)cf>(y)Sp(cT)Vs(x, y, tr) 
| J d D xd D ydcrd~i ^(x)^(y)5p( < T)5p( 7 )V4(x, y, a, 7 ), 



(113) 



where 



y 3 (x, y, <t) = [(z - p/(0))*(x - y)8{a - x) - p<5(x - y)/(x - «r) - p/(x - y) (<5(x - <r) + 5(a - y)) 
^4(x,y,<r,7) = pS(j-x)f(x-y)S(y - tr) - p£(7 - x)<5(x - y)/(x - tr). 



(114) 
(115) 



Note that the first term of Eq. (|112[) . when computed in the limit Sp = of the action (Eq. (|113[1 ) yields the bare 
propagator G (p, z). The fist term then corresponds to the free (gaussian) part of the field theory and the terms 
involving three and four fields are the interacting part. The latter can be treated perturbatively using standard 
diagrammatic techniques of field theory. One can easily see from the form of the interacting terms that in such 
diagrams no loops involving the dp field may arise because, at fixed disorder, it acts as an external field while a 
generic n-loop diagram comes from the average over the disorder and yields a l/p n contribution to the resolvent. 



A. The correlation functions for the density field 

In order to perform the loop expansion one needs the n-point correlation functions of the Sp field. Since 

1 N 



^=i + w=/(nV^£ d >- 



Xk) = - P v =1 ' 



(116) 



one has 



5p(x) = 0. 



(117) 
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Similarly, the fact that 



p(x)p(y) = 1 + 5p(x) + 5p(y) + Sp{x)5p(y) = 1 + 6p(x)6p(y) = 

N N „/ N , D 



i r( N h^x A i N r ( N d D x \ (118) 

7 E /(n ^ J * (x - Xfe)5(y -^vE/(n^j '(* - **)*(y - **) 



JV(JV-1) N . 1 1 , 



implies that 



Sp(x)<5p(y) = -<5(x-y). (119) 



To carry our the perturbative expansion up to order 1/p 2 , the 3- and 4-point correlations are needed. These can be 
derived according to the lines described above, giving 

6p(x)8 P (y)6p(z) = 4^(x - y)<5(y - z) = 6p(x)5p(y) ■ 6p(y)Sp(z), (120) 
P 

5p(x)6p(y)6p(z)6p(t) = -I J(x - y)6(y - z)«5(z - 1) + 6p(x)6p(y) ■ 6p(z)6p(t)+ (121) 
p 6 



S P (x)5p(z) ■ Sp(y)Sp(t) + 5p(x)6p(t) ■ 8p(y)Sp(z). 



1. The general expression 



We may write as well the expression for the arbitrary n-point correlation Sp(y 1 )Sp(y 2 ) ■ ■ ■ Sp(y k ), needed to compute 
the self-energy to order 1/p 3 or higher in the field theory. To give our result, we shall need some notations. 

Let u> be an arbitrary partition of the set {1,2, ... ,k} into subsets. For instance, for k — 4, uj could be a par- 
tition into two subsets, such as w = { {1, 2}, {3, 4}} or u — { {1, 3}, {2, 4}}, or a partition into 4 subsets such as 
{ {1}, {2}, {3}, {4} }, etc. Let ||a|| be the cardinality of the set a, for instance, if to = { {1, 2}, {3, 4}}, then ||o;|| = 2 . 

We also define T>( k \ the set of all possible partitions of {1,2,..., k}. Given a partition w, the subsets associated to 
it will be J?;, w , with I = 1, 2, . . . , ||w|| . We shall need to consider %( k >, a subset of the set of all partitions V^. 
is made of all partitions u> such that > 1 for all Z = 1, 2, ... , ||cj||, i.e. partitions in which none of the subsets 

contains less than two integers. Then the general result is: 

Sp(yi)Sp(y 2 )...Sp(y k ) = 

The proof is given in app. To recover Eq. (|120p from this formula, note that the set W 3 ' of allowed partitions of 
{x, y, z} contains a single partition, with just one subset (||w|| = 1), namely u> = {{x, y,z}}. On the other hand, to 
obtain Eq. (|12ip we need the set of allowed partitions for {{x, y, z,t}}. There are four such partitions, namely 
wi = {x,y,z,t}, uj 2 = { {x,y},{z,t} }, uj 3 = { {x, z}, {y, t} }, and uj 4 = { {x, t}, {y, z} }. Clearly, ||wi|| = 1, while 

B. Diagrammatic expansion: one loop 

In order to write down the one-loop term it turns out to be convenient to write V3 and V4 in terms of the interaction 
vertex (Eq. (|Tg|) ): 

jd D xd D a^ +i ^V 3 (x,y,a) - [G^( Pl ) - V(p 2 , - Pl )] e^+^> y = /i( Pl , p 2 )e^ +p ^ (123) 
y'd I5 xd I) crd I) 7e ip i x+i ^ <T+ip 3^y 4 (x,y, ( T,7) = -y(p 2 , - Pl )e ?;(Pl+P2+P3)y . (124) 





(122) 
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-G (p,2) 
P 




FIG. 1: Diagrammatic notation 

The latter expression depends only on two momenta. Tims when the vertex V4 is involved, one has to make its 
expression symmetric by joining the Sp propagators with the two possible external links offered by this vertex. 

Fig. Q] defines our diagrammatic notation. Note that the vertex V4 is not symmetric under the interchange of the 
two bp lines. Now we are able to write down the one-loop diagrams. 

Recalling that the resolvent G(p, z) is given by Eq. (|112p . we compute the one-loop contribution to ((/>(x)(/>(y)): 



4 1J - /^o(p-q,^ 2 (p,-q) (125) 



= (126) 



The last diagram gives a general result: every tadpole made with a vertex with four fields gives a vanishing contribution 
due to the form of the vertex. The term with one external Sp insertion, arising from <5p(x) (<^>(x)</>(y)), is given by 



r (D 2G 2 (p,z) f d D q 

L\ = c • = / 7^-TdG (p- q,z)^(p,-q) (127) 

p J {2n) 1J 



The last contribution to the self-energy at one loop comes from 5p(x.)5p(y) (c/>(x)(/>(y)), and is 




4 1} ^ 2 = \ J (0T G o( P - q, z) (128) 

Note that this last contribution has an ultraviolet divergence since the propagator goes to a finite costant when the 
internal momentum goes to infinity. Nevertheless, by adding the four diagrams 

L« + L« + 4 3) = - J (0jG o (p - q, z) [/i(p, -q)G (p, z) - if = (129) 

(130) 



J ^Go(d,z)V 2 ( (l ,p)^Gl( P ,z)S^(p,z), 
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the divergence disappears and one recovers the combinatorial result for the one-loop self-energy. 

C. Two loops 



Let us first consider the two-loop diagrams arising from (</>(x)</>(y)) : 



(2)_ £ I _ Gl(p,z) f d^qd^k , 



LY> = 4 » » » = / " gg(q, z)G (q - k, z)^(p, q - p)^(q, -k) (131) 



r(2) 
^2 



/^T^i = 9MPlll f d ^ q ^ k G (q, z)G (k, z)G (p q k, z) M (p, q - p)/i(p, k - p). (132) 

• /"(q>k- p)/x(k,q- p) (133) 

if = = _^i£) | d ( ^p k Go(q, z)G (k, z) M (p, q - P)M(P, k - q)MP, k - p) (134) 

/^" - T^i = | d ( ^p" Go(q, s)Go(k, z)/i(p, q - p) M (P, k - p). (135) 



r(2) 
^4 



[F(p-q,-q)-F(k-p,-q)] (136) 



4^ 4* 1 = _gg^f) | i_S^ Go(q , z)Go(q _ kjZ ) Mq ,_ k ). (137) 

■[-y(p-q,k-q)-y(k,k-q)] (138) 




4 2) = — * -4 — = / d ( 2^2D k G ° 



(139) 



^ EE f \ f\ = I -^Sp Go ( k , z ) Go ( q , z)A1 2 (piq _ p)Ai2(pjk _ p) (14Q) 

(2) 

The diagram L 7 seems to be already included in the Dyson resummation of the one-loop result. However, since 
diagrams with one and zero external legs have to be included in the diagrammatic expansion, it also provides a genuine 

C2) _____ 

contribution to the two- loop result. Note that in order to obatin L ( 3 ' one uses Eq. (|120[) for the 3-point correlation 
of bp. The other diagrams involve only the disconnected part of the 4-point correlation function, while the connected 
one would only matter at three loops. 

Next we must copsMervthc contribution arising from <5p(x) (0(x)0(y)): 

2 Go (p,z) [d D (id D k 




P 2 



f d^qd^k 

J ( 27r )2-D G o(q. z)G (q - k, z)p(p, q - p)^ 2 (q, -k) (141) 



4^ = < ^^\ = - 2G0{ J> ,Z) J ^^G (q, z) Go (k, z) Go (p q k, z)p 2 (k, q - p> (142) 



/Lt(q,k-p)/x(p,k-p) 



(143) 
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L$^ £—£^i = ggg(Pif) | ^Sp Go(q , z)Go(kjZ ) Mq , k _ q ) M p, k _ p) (144) 



Li? = ^ i ^ = - 2G0{ *' Z) I d ^2D G ' 



(145) 



r(2) 
^12 



fJ ^V^4 = - 2G °^ 2 P,Z) | d ( ^p* Go(q, «)Gf (k, Z>(p, k p) [V(p - q, -q) + V(k - p, -q)] 



(146) 



L^ ^W^ =-g^Lg) / ^pG ( q ,z)Go(k I z) / , 2 ( P ,k- P )Mp,q-p) (147) 
As before, we have used the disconnected part of the 4-point function, appart from L 10 where we have used the 

(2) 

3-point function. Note also that L 13 arises both in the Dyson resummation and in the two-loop expansion. 



Finally, we consider the diagrams arising from 6p(x)5p(y) (</>(x)</>(y)): 




L$ = <- i = 1 1 ^^ G 2( q ,z)Go(q-k I z) / , 2 (q,-k) (148) 

L< 2) = f^^ i ^ d ( ^^ k G (q, z)G (K z)G (p - q k, z) M 2 (q, k - p)/i(k, q - p) (149) 

L$ = ^^V^ = "4 / d ^ d ^G (q,z)Go(k,z) Ai (q,k-q) (150) 

L<$ = = ± J d ^ d ^G (q, z)G„(k, z) [V{p - q, -q) + V(k - p, -q)] (151) 

L { ^= ^\ ^ =j 2 J ^^ Go (q, z)G (K z)n(p, q - p)M(p, k - p) (152) 

We now show how the diagrams can be summed up to give the combinatorial expressions for the self energy. 
Consider the diagram L^; it has the same topology (in the sense of momenta flow and vertex positions) of . In 

(2) (2) 

fact, it can be combined with L y s ' and L\l to give: 

L< 2) + 4 2) + L< 2) = 1 j d ( ^ d ^ k G 2 (p, z)G 2 (q, z)Go(q - k, z)V 2 (p - q, p) M 2 (q, -k) = (153) 
= 4 2) (p,z) + l2 1 (p,z), (154) 

where we have defined 

1 f d^ad^k 

«i(p, z) = J * [G 2 (p, z)G (k, z)y 2 (p q, p) 2G 2 (p, z)G (q, z)G Q (k, z)V 2 (p - q, p)F(q - k, qj\ . 

(155) 



23 



In the same way we can combine L\ , h\' and L\ 5 , they have the same topology of Sb(p, z): 

Lf + + = j 2 J d ( ^^ k G (q, z)G (k, z)G (p - q - k, z) M (q, k - p) M (k, q - p)- (156) 

• V(p - q, pMp - k,p) = 4 2) (p, z) + fl 2 (p, «) (157) 

where 

1 f d^ad^k 

«a(p, *) = ^ _/ r^p [ G o(P, ^)G (P - q - k, z)V(p - q, p)V(p - k, p)+ (158) 

-2G§(p, z)G (q, ^)G (k, z)V(q + k, p)f (p - k, p)V(q + k, k)] . (159) 

(2) (2) (2) 

We now add the diagrams L\ , L\ 3 ' and L^g because they produce the Dyson resummation of the self energy at 
one loop that we want to isolate from the other contributions that have to be included in the self energy at two loops. 
They give 

4 2) + 4 2) + Lg> = i J d ( 2^2D k Go(P> z)G (q, z)G (k, z)F(p - q, p)V(p - k, p)/i(p, q- p) (160) 
■ /x(p, k - p) = G (p, z)^ (1) (p, z)G (p, z)S (1 \p, z)G (p, z) + f2 3 (p, 2), (161) 

where 

1 f d^ad^k 

Os(p, Z ) = - 2 J * [G (p, z)G (q, z)G (k, z)V(p - q, p)V(p - k, p) (1 - 2G ( P , z)V(p - q, p))] . (162) 

At this point one can check that 

Q 1 (p, z) + I2 2 (p, z) + fl 3 (p, z) + Lf + hf + 4 2) + 4 2) + L$ + + Z$ + Z$ + L<$ = £™ (p, z) (163) 
and the combinatorial result is recovered. 



D. The small p behaviour 



We will now prove that the prefactor of the term p 2 \( D ~ 2 y 2 is zero to all orders in perturbation theory. For this 
purpose, the field theory approach turns out to be very convenient. 
Consider the vertex V3 with three fields, 

= Go \p, z) + V(q, p) = z - pf(0) + pf(p) + pf(q) - pf(p - q). (164) 

It is easy to see that this vertex is symmetric and can be written as 

V 3 = z + S(p,q), where S(p, q) = 5(q, p). (165) 
Moreover one can check directly that 

5(p, 0) = 5(0, p) = 0. (166) 

Consider now the vertex with four fields. We see that the Wick contractions between the fields 8p symmetrize the 
vertex. In fact in every diagram this vertex appears in the form 

= -(V(k,-p) + V(q-p-k,-p)) (167) 
= -p (/» - f(P + k) + f(q - P - k) f(q - k)) . (168) 
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The important thing is that this vertex vanishes when one of the two Go bare propagators carries a null momentum. 
Consider now a diagram that arises from the expansion of the resolvent G(p,z). At the lowest order in z when the 
diagram contains some three-field vertices one has to consider only the symmetric part of these vertices. Let us 
apply the method explained above in order to extract the contribution to the self energy proportional to z( D ~ 2 ^ 2 . 
Apparently, if one sets to zero the momentum of a bare propagator that enters into a vertex then its contribution to 
the imaginary part vanishes. This seems very strange because from this argument it follows that only Lg 1 ' contributes 
to the imaginary part. Moreover if we consider the two-loop contributions we see that there are no contributions 

(2) (2) 

to the imaginary part of the self energy because the diagrammatic expansion L\ -L\ s contains at least one vertex 
that vanishes when we set to zero one of the momentum brought by a (/)-propagator. However the argument is not 

(21 (2) 

complete. Actually, the diagrammatic expansion L\ -L\% does not contain only the contribution 

G (p,z)U^(p,z)G (p 7 z) (169) 

since it contains also the Dyson rcsummation of the one loop self energy. This is the fact that completes the argument 
and will lead us to prove that a contribution proportional to 

z (D-2)/2 p 2 

cannot appear at any order in perturbation 

theory. 

We start checking the argument just given at the one-loop level. Let us introduce the notation 

=G - 1 (p,z), (170) 

so that the self energy at one loop can be written diagrammatically as 



S^(p,z) = k •+ • 1 + (171 



From this expansion and from the above argument one sees that the imaginary part of the self energy (we will refer 
always to the imaginary part proportional to ^( B ~ 2 )/ 2 ) may come from the last diagram and is correctly given by 

Im lim E (1) (p, z + ie) oc z (D - 2)/2 V 2 (p, p), (172) 

e->0+ 

which is also the contribution that can be easily calculated from the combinatorial expression. Now consider the 
expansion at two loops. From the combinatorial expressions of the self energy we immediately see that the immaginary 

(2) 

part comes from only S c and can be rewritten in the form 

Imlim£ (2) (p,z + ie) oc -2z (D - 2)/2 V(p, p)S (1) (p, z). (173) 

(2) (2) 

Consider now the diagrammatic expansion for the two loop self-energy L\ -L\ s . We have to extract from this 
expansion the term 

S^(p,z)G {p,z)S^{p,z). (174) 

Now we will do this in a diagrammatic way. Consider the diagrammatic expansion for the above term (Fig. [5]). 
If we want the imaginary part of the self energy proportional to z^ D ~ 2 '' 2 we have to consider the term 

18 

A(p, z) = [G (p, z)}- 1 L? ] [G (P, z)]~ l 

i=l 

and the diagrams in Fig. O When we calculate this contribution we have to set to zero the momentum carried by 
one internal propagator Go so that the contribution coming from A(p, z) does not matter. We have to calculate only 
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CD, 



+ 



+ 



(175) 
(176) 
(177) 

(178) 
(179) 



FIG. 2: Diagrammatic expression for Eq. (|174[) 



the term coming from the Dyson resummation so that the imaginary part of the self energy at two loops is given by 



Im lim £ (2) (p, z + ie) = -2 



Im lim E^\p, z + ie) 



2i7 (1) (p,z)Im lim 

6->0 + 



Im lim 



-2I7 (1) (p,z)Im lim 

e-»0+ 



OC 



_ 2z (D-2)/2 y(pjp)i; (l) (p5z) 



At this point we can give also the analytical argument 

Im lim £" (2) (p, z + ie) = 

e— >0 

= -Im lim £ (1) (p, z + ie)G (p i z + ie)Z' (1) (p, z + ie) 
<x-2V(p,p)ZW(p,z) 

where we have used the fact that 

E(0, z) = . 



(180) 
(181) 

(182) 

(183) 

(184) 
(185) 
(186) 

(187) 



On the same line we can give the imaginary part proportional to z^ D 2 ^ 2 at three loops because this contribution 
comes from the Dyson resummation of one and two loops self energy: 



ImlimX'( 3 )(p ! 2 + ie) = -Imlim Z^G E (1) Gq£ (1) + 2S^G E^ 

e— >0 e—>0 . 



oc z^' 2 3 



2V(p,p)Z ( V 



(188) 
(189) 



At this point we can give a general expression for the imaginary part proportional to z^ D 2 ^ 2 at any perturbative 
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order: 



Imlimi: (Tt) (p,z + ie) = 



e-s-0 



-Im lim 

e->0 



n k 

E E II (^ (ia) (p^ + ")Go(p^ + ic) 



k—2 a— 1 



[GoCp^ + ie)]" 1 ] 



(190) 

(191) 
(192) 



From this expression we can prove by induction that the imaginary part of the self energy proportional to z^ D ^ 2 ^ 2 
cannot appear at any order in perturbation theory. In fact we have seen that it docs not appear at one and two loop 
so we can prove that if it does not appear up to n loops it does not appear ton+1 loops too. We can see that 



Im lim E (k) (p, z + ie)G (p, z + ie) oc z (D ~ 2)/ V 
where 7 > 2 and where we have showed only the term at lowest order in z. Moreover we have 

£W(p,z)G (p,z)~l + O(p 2 ) 
where we have neglected the higher order in z. It follows that the generic term in (|192p is of order 



(193) 

(194) 
(195) 



with (3 > 4 because [Go(p, ; 



p 2 . This completes the proof. 



VI. CONCLUSIONS 



In conclusion, we have given a detailed description of the perturbative high-density computation of the resolvent (and 
in particular the density of states) of Euclidean random matrices within two different formalisms. The combinatorial 
formalism of scc. HVI results in fewer diagrams and is probably more convenient when the goal is to obtain an expression 
of the self-energy at a given order. On the other hand, the field-theoretic formalism (sec. [Vj) , though producing a 
higher number of diagrams, has allowed us to analyze the p — > behavior at all orders in perturbation theory. This 
analysis shows that the immaginary part of the self-energy in the limit of small momenta (which controls the width 
of the Brillouin peak of the dynamic structure factor) has, in contrast to previous claims [l5l - [l8l . l30j . the structure 

2 

-ImX!(\,p) =B\ E * 1 p 4 +C\%^ + ... , (196) 

cr 

where C,B > are amplitudes, and c is the speed of sound. This implies in particular a p 4 scaling for the Brillouin 
peak width, but it also shows that the structure of the theory is more complex than in the case of scattering from 
lattice models (29|. 
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Appendix A: Proof of Eq. (Tl22l) 



The proof proceeds by induction. First note that the explicit computation in Eqs. (|117[) and (|119[) already implies 
that Eq. [TH] holds for k = 1 and k = 2. 

The cornerstone of the proof is a general result for the fc-point correlation functions of p (rather than Sp). The 
sought correlation function, in the thermodynamic limit, is 



Piyjpfa) ■ ■ ■ P(Yk) = !+ E 



1 



fc-IMI 



i>€V<- k \\\bj\\<k 



P 



n 

j=i,h«j,ji>i 



(Al) 



Eq. (|A1I) looks very similar to Eq. (|122[) . yet we note the following crucial differences: 
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• The partitions us belong to V^ k ' rather than to the restricted set T-L^ k \ In particular, the term equal to 1 in 
Eq. (|Aip follows from the only partition us with |jw|| = k, namely { {1}, {3}, . . . , {k} }, which obviously does not 
belong to % {k \ 

• In the innermost product in Eq. (|A1[) . a subset fi^ u with just one element, ||^?; itlJ || = 1, merely contributes a 
factor of one. Hence, for all practical purposes, such a subset £li^ can be ignored. 

To establish Eq. (|Al[) . first note that 

1 N 

p(yi)p(y 2 )---p(y/e) = — Yl ^ x n)^(yi - x i 2 ) - x *Ji ( A2 ) 

where the average is taken with respect to the fiat probability measure, 

Now, for a given assignment of the k particle labels i\, £2,. ■ ■ ,ik, we declare that all terms with a coinciding particle 
label i r form a subset It is then obvious that every assignment of the k particle labels i\, «2,- ■ • 4k defines a 

partition us in V^ k \ Furthermore, a little reflection shows that all possible partitions in V^ k ' can be obtained in this 
way. Eq. (|A1[) follows from the following three facts about a given partition us: 

1. There arc N(N — 1) . . . (N — possible assignments of the k particle labels i\, t2,. ■ ■ ,ik that yield the 
partition us (you are given N choices for the particle that appears in the subset N — 1 for that appearing 
in i?2,u, and so forth). 

2. A subset with a single clement, = 1, contributes a factor l/V. 

3. A subset with more than one element, ||^z,w|| > 1, contributes a factor 

r— 1 



Now consider a partition us that belongs to but does not belong to . Imagine that us contains k — s subsets 
£2i lU with just one element. The values that s can take are s = 0, 2, 3, 4, . . . , k — 1. We are not interested in the trivial 
case s = 0, that corresponds to the partition { {1}, {3}, . . . , {k} }. Hence, for s > 0, we simply erase from us all the 
k — s subsets f2i_ u with ||^ jS || = 1. The s integers 

A = {ptPt---,Ps}, 

that belong to the remaining ^2;.^, form the irreducible set A associated to the partition us. The list of the fli iU1 with 
11^,^11 > 1 provides a partition us of A, that obviously belongs to H( S ' A K Furthermore, \\us\\ = k — s + |cj||, so we have 
for the prefactor in Eq. (|A1I) that 

1 f 

pk-\\ui\\ ~ pS-\\&\\ ' 

Hence, since s < k, the induction hypothesis implies that the added contribution in Eq. (|A1[) of all the partitions 
sharing the same irreducible set, A, is 

At this point, we may rewrite Eq. (|A1[) as 

p(yiMy 2 )---p(yfc) = i+ Y -^p, 

fe-i 

Y E ^(y^ ) <5 p(y^ ) ■ • ■ § p(y^ ) ■ 



INI 

n n ^ 
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We finally note that, if one writes p(y r ) — 1 + Sp(y r ), 



P(yi)p(ya) • • • P(yfe) = 1 + Sp(y!)5p(y 2 ) . . . Sp(y k )+ 



fe_1 (A4) 

51 s p(yp^) 6 p(y^)--- s p(y^)- 



s=2 A={P*,l3f,...,pA} 

Comparison of Eqs. (|A3|I and (|A4|) completes the proof. 
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